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Using data taken with the CLEO II detector, we present the first full angular 
analysis in the color-suppressed modes B^ —^ J/tp K*^ and B~^ — > J/ifjK*^. 



C^ ' This leads to a complete determination of the decay amplitudes of these modes. 



In addition, we update the branching fractions for B — > J/ip K and B 
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One of the interests in B —>■ J/ip K* decays is their role in CP violation measurements 
at asymmetric -B-factories. The vector-vector decay B^ -^ J/ipK*^, with K*'^ -^ Kgir^, is 
a mixture of CP-even and CP-odd eigenstates since it can proceed via an S, P, or D wave 
decay. If one CP eigenstate dominates or if the two CP eigenstates can be separated, this 
decay can be used to measure the angle f3 of the unitarity triangle in a manner similar to 
which the CP-odd eigenstate 5° -^ J/ip K^ is used. 

Measurements of the decay amplitudes of 5 — *> J/i}) K^*^ transitions also provide a test of 
the factorization hypothesis in decays with internal VT-emission. Several phenomenological 
models, based on the factorization hypothesis, predict the longitudinal polarization fraction 
in B —>■ J/tp K*, denoted Fl/F, and the ratio of vector to pseudoscalar meson production, 
R = B{B -^ J/^K*)/B{B -^ J/ijjK) §^. It has been noted [|,|, that form factor models 
cannot simultaneously explain the earlier experimental data for these two quantities. The 
high values of Fl/F measured by ARGUS and CLEO II |§|, with low statistics, are 
not consistent with factorization and the measured value of R. The CDF collaboration has 
measured a lower value of Fl/F 0. Additional information about the validity of factorization 
can be obtained by a measurement of the decay amplitude phases, since any non-trivial phase 
differences indicate final state interactions and the breakdown of factorization [p!0| . 

In this paper we present a complete angular analysis and an update of the branching 
fractions for B — > J/^ K^*'^ decays using the full CLEO II data sample. Assuming isospin 
symmetry, we determine the fraction of longitudinal polarization, the parity content and 
the phase differences of the decay amplitudes from the modes B^ -^ J/ip K*^ and B^ — > 
J/ipK*^ using the K*+ and K*^ decay modes to K+tt^, K^7c+, K+tt", and i\:°7r°. The J/tp is 
reconstructed in its leptonic decay modes to e~^e~ and n'^n~ ■ The measurements presented 
here supersede previous CLEO II results |§] , which are based on a subset of the data used 
for this analysis. 

The decay B -^ J/ipK* is described by three complex decay amplitudes. Following a 



suggestion of Dunietz et al. [|1T],|I2|, we measure the decay amplitudes Aq = —Jl/3S + 



2/3 D, All = V 2/3 S + \/l/3 D, and A^_ = P, where S, P, and D denote S, P, and D wave 

II 



amplitudes, respectively. Normalizing the decay amplitudes to |AoP + l^llP + I^J-P = 1 and 



eliminating one overall phase leaves four independent parameters. 

The full angular distribution of a -B meson decaying into two vector particles is specified 
by three angles. Previously the helicity angle basis [|l3l has been used for angular analyses 



oi B ^ J/ip K* decays. Because of its convenience for extracting the parity information, we 
use a different set of angles, called the transversity basis |T^ . The direction of the K* in the 
J/ip rest frame defines the x-axis of a right-handed coordinate system. The Ktt plane fixes 
the y-axis with Py{K) > and the normal to this plane defines the z-axis. The transversity 
angles 6'tr and 0tr are then defined as polar and azimuth angles of the l~^ in the J/ip rest 
frame. The third angle, the K* decay angle 6k*, is defined as that of the K in the K* rest 
frame relative to the negative of the J/ip direction in that frame. Using these definitions the 
full angular distribution of the B -^ J/ip K* decay is |]12| : 






1 d^r 

r d cos 6'tr d cos 6k* d^tr 

= ^{2\Ao\^ cos2^i^.(l-sm2etrCos2 
+ |y4|||^ sin^6'K*(l -sm^6'trsin^ 
+ |v4_l|^ sin^ 9k* sin^ 6'ti. sin^ 0tr 
— Im (yl|Jylx) sin^ 9k* sin 2^ti. sin 0tr 
+ ^ Re (Ap A|| ) sin 261;^* sin^ 9tr sin 20tr 
+ 4= Im {AqA±) sin 26'^^. sin 26'tr cos 0tr }• 

For i? decays the interference terms containing A± switch sign while all other terms remain 
unchanged. 

The data for this analysis were recorded with the CLEO II detector located at the Cornell 
Electron Storage Ring (CESR). We have used a data sample of approximately 3.4 x 10^ BB 
events taken on the T{4S) resonance and representing an integrated luminosity of 3.1 fb~^. 
To evaluate non-66 backgrounds, we have also collected a "continuum" data sample 60 MeV 
below the T(4S) resonance, with an integrated luminosity of about 1.6 fb~^. 

The components of the CLEO II detector [Q most relevant to this analysis are the 
charged particle tracking, the Csl electromagnetic calorimeter and the muon counters. The 
tracking system comprises a set of precision drift chambers totaling 67 layers inside a 1.5 T 
solenoidal magnet. It measures both momentum and specific ionization {dE/dx) of charged 
particles. 

Electron candidates are identified by their energy deposition in the calorimeter, which 
must equal their measured momenta, and their specific ionization, which must be consistent 
with that expected for electrons. At least one muon candidate is required to have penetrated 
five nuclear interaction lengths of material while the other must have penetrated at least three 
interaction lengths. The decays B^ —>■ J/ip K^ and B^ —>■ J/ip K^ have little background, 
therefore only one of the two leptons has to be positively identified. In this case the identified 
muon candidate must penetrate at least three interaction lengths. We require the dimuon 
invariant mass to be within 45 MeV/c^ of the J/ip mass, which corresponds to a 3o" selection. 
For the dielectron invariant mass we require —150 MeV/c^ < mee — trij/^ < 45 MeV/c^ to 
allow for the radiative tail. The J/ip energy resolution is improved by a factor of 5 — 6 by 
performing a kinematic fit of the dilepton mass to the nominal J/ip mass. The kinematic 
fit does not affect the measurements of the transversity angles. The resolution of the angle 
measurements is better than 0.06 radian for all decay angles. 

We require the charged hadron candidates to have dE/dx measurements that lie within 
three standard deviations (a) of the expected values. We reconstruct Kg candidates through 
the decay to tt+tt^ and we reconstruct 7r° candidates through the decay to 77. Candidate 
K* mesons are required to have a Ktt invariant mass within 75 MeV/c^ of the nominal K* 
mass. 

In symmetric e+e^ annihilations at the T(4S) resonance, the energy of a -B meson must 
equal the beam energy. We require the energy difference \AE\ between the B candidate 
and the beam energy to be less than 45 MeV for J/ip K^ and J/ip K^, less than 30 MeV for 
J/^ (K+n-) and J/^ (i^s7r+), and less than 60 MeV for J/ip {K+n°) and J/ip (i^^vrO). These 



ranges correspond to approximately 3o" in jAE"!. In the rare case, that an event has more 
than one candidate per mode we keep only the candidate with the highest probability, which 
is based on the measured \AE\ and, if available, the measured dE/dx, the 7r° mass, and 
the time-of-flight information. The resolution on the beam energy is an order of magnitude 
better than the resolution on the B candidate energy. Therefore we substitute the beam 
energy in the calculation of the -B-candidate mass (referred to as the "beam- constrained 
mass" tub)- The detection efficiencies range from 48% for the B^ — > J/ipK^ mode down to 
9% for 5° -^ Jji) K*° with K*° -^ K^ti^. 

The most severe background in the B -^ J/ipK* modes is feed-across from one B -^ 
J/il)K* mode to another. For such events both the total energy and the beam-constrained 
mass are very close to the signal region. The biggest source of feed-across background is 
from swapping a random or misidentified slow vr^ for the correct one. Consequently most 
background events have the vr" moving backwards with respect to the K* direction of flight. 
To suppress this background we require the K* decay angle to satisfy cosOk* < 0.7 in 
these decays. This is equivalent to a constraint on the vr^ momentum, corresponding to a 
minimum p^o of about 200 MeV/c. The total fraction of feed-across events in the signal 
region, averaged over all K* modes, is 8.0%. 

There might be a contribution from non-resonant B -^ J/ip Ktc decays in the K* signal 
region, though neither the previous CLEO measurement ^ nor CDF found events in the 
K* mass sidebands. However, examining the Ktt invariant mass spectrum (Fig. [^) shows an 
excess of events between 1.1 and 1.45 GeV/c^. By computing the kinematics of non-resonant 
B -^ J/ipXs decays, using both the J/tp momentum spectrum from inclusive B decays |jl5| 



and several theoretical models |T6[, we do expect strangeness-containing final states with 
invariant masses in this region. Decays via higher K* resonances may have line shapes 



consistent with the rriKn distribution seen by us ^7j. Unfortunately, due to the limited 
statistics for rnxn > 1-1 GeV/c^, we cannot distinguish between possible components. In 
addition, we considered many other possible origins for these events including feed-down 
from modes such as -B — > J/ip K*7t or i? — >■ J/ip Kp, and feed-through from B -^ J/ip K, 
and found none of these to contribute significantly. We estimate the amount of the non- 
K*{892) contribution in the signal region to be 6.4% with a conservatively chosen systematic 
uncertainty of ±100%. 

Other misidentified B -^ J/ip X decays, like B —>■ J/ip K*7r, do not contribute significantly 
to the background since they lie outside the energy window. With a similar analysis CLEO 



has found 9 events for B^ — >■ J/ip p° [ffH] . If a pion from the p° is misidentified as a kaon, 
rriKiT could fall in the K* region but these events would fail the \AE\ energy criterion. Feed- 
across between the J/ip K and J/ip K* modes is also suppressed by the requirement on | AE'l. 
Furthermore, the contributions are uniform in the beam-constrained mass. 

We define combinatorial backgrounds to be events that do not contain a true J/ip -^ l^l~ 
decay. In both the BB Monte Carlo simulation and our continuum data sample we see very 
few such events. 

We must correct our data for detection efficiency. To obtain the efficiency as a function 
of all three angles, a large Monte Carlo sample (120,000 events/i^* mode) is divided into a 
20x20x10 grid in cos6'tr, cosOk* and 0tr. For each J/ip K* final state the efficiency is fitted 
separately with polynomials in three dimensions including all correlations. The efficiency 
distributions are nearly uniform in all angles except the K* decay angle, where it drops at 
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FIG. 1. The niKn distribution for niB > 5.27 GeV/c^. The data points, the fitted K*{892) 
mass peak including feed-across (empty histogram), and the combinatorial background (shaded 
histogram) are shown. 



high cos 9k* because of the slow pion. 

To determine the decay ampUtudes, a four-dimensional unbinned maximum likelihood fit 
is performed to the distributions of the three angles and the beam-constrained mass. Setting 
(p{Ao) = 0, we fit for the longitudinal polarization fraction, |v4op = Tl/F, the parity-odd 
fraction, |v4_Lp = |Pp, and the phases (f){A\\) and (f){A±). Other free parameters in the 
fit are the branching fraction B{B —>■ J/tpK*), the mean of the itlb distribution and the 
normalization of the combinatorial background of each mode. Fitting for the branching 
fraction and the polarization parameters simultaneously ensures the correct treatment of the 
background events and automatically adjusts the branching fraction measurement for the 
polarization dependence of the efficiency. The one- dimensional projections of the resulting fit 
function are shown in Fig. ^j. The results are listed in Table |. The correlations between the 
fit parameters are small. The systematic uncertainties of the decay amplitude measurements 
are dominated by those in the efficiency parameterization and background polarization and 
are small compared to the statistical errors. 

We also repeated the fit to the decay amplitudes using helicity angles rather than transver- 
sity angles as well as performing one-dimensional fits to both the longitudinal polarization 
fraction and the parity-odd component. An independent angular analysis with the same 

PI 



to evaluate the 



data sample has also been performed, using a Monte Carlo technique 
likelihood function. All results are in agreement with those reported here. 

These results are the first determination of the parity-odd component and the phases of 
the decay amplitudes of the B -^ J/ip K* decay. The small fraction of the parity-odd compo- 
nent encourages using the B^ -^ J/ip Kgir^ decay for CP violation studies at asymmetric B 
factories. The phases of the decay amplitudes are measured to be close to zero or vr, giving 
no evidence for strong final state interactions. 

The branching fractions of the B -^ J/ipK* decays are a result of the angular fit. To 
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Parameter 



Value 






0.52 ± 0.07 ± 0.04 
0.16 ± 0.08 ± 0.04 
-0.11 ± 0.46 ± 0.03 radian 
3.00 ± 0.37 ± 0.04 radian 



TABLE I. Resulting decay amplitudes from the fit to the transversity angles. The phase 0(^o) 
has been set to zero. The first error is statistical and the second is the estimated systematic 
uncertainty. 
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FIG. 2. One-dimensional projections of the four-dimensional fit to the B -^ J/tp K* data. The 
plot of the beam-constrained mass shows the data (histogram), the fit (solid line), the sum of all 
backgrounds (dashed) and the B -^ J/tp K* feed-across (dotted) . The angular distributions are 
background subtracted and efficiency corrected. 



measure the B —>■ J/ipK^ /K^ branching fractions we performed one-dimensional fits to the 
beam-constrained mass distributions. The results of the fits are shown in Fig. ^. All measured 
branching fractions are listed in Table |I|, where we have assumed that the production rate 
of neutral and charged B mesons is the same on the T(4S') resonance, in agreement with the 
actual measured value of /+_//oo = 1.12 ± 0.20 [^ and a theoretical prediction ^T|]. The 
main sources of systematic uncertainties of the B -^ J/ip K branching fraction measurements 
are track finding, track fitting, and lepton identification efficiencies, and the uncertainty of 
the world average of B{J/ip — > /"'"/") [^. In the B —>■ J/ip K* branching ratios, uncertainties 
in the amount of feed-across and non-i^* decays dominate the systematic error. 
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Decay mode 



Signal Yield 



Branching Fraction [10 



B+ -^ J/il) K+ 

B^ -^ J/i) K^ 

B+ -^ J/i) K*+ 

B^ -^ J/i) K*^ 



198.1 ± 14.9 



45.5 



+7.3 



-6.6 

42.5 ± 7.1 
81.6 ± 10.3 



1.02 ± 0.08 ± 0.07 

0.85 +°;^2 =t 0-06 
1.41 ± 0.23 ± 0.24 
1.32 ± 0.17 ± 0.17 



TABLE II. Measured signal yields and branching fractions. 



With the assumption of equal partial widths, T{B+ -^ J/^ K'^*^+) = T{B° -^ J/ip K^*^°) 
and eliminating common systematic uncertainties we determine 



/+- ^B+ 
/OO TfiO 



1.15 ± 0.17 ± 0.06, 



Assuming isospin invariance, we find for the ratio of pseudoscalar to vector meson pro- 
duction 



B{B ^ JibK*) 

R = -1 '-L 1 = 1.45 ± 0.20 ± 0.17. 

B{B -^ J/i) K) 

These new measurements of Fl/F and the ratio R indicate that the discrepancy with 
naive factorization models is not as acute as before. 
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